QED Calculation of E1M1 and E1E2 Transition Probabilities in One-Electron Ions 
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The quantum electrodynamical (QED) theory of the two-photon transitions in hydrogenlike ions 

■ is presented. The emission probability for 2s 1 / 2 — * 27(_E1) + ls^ transitions is calculated and 
compared to the results of the previous calculations. The emission probabilities 2pi/2 — * j(El) + 
7(-E2) + lsj/2 and 2pi/ 2 — * l(El) + j(Mi) + ls^ are also calculated for the nuclear charge Z 
values 1 < Z < 100. This is the hrst calculation of the two latter probabilities. The results are 
given in two different gauges. 

<: 

1. Introduction. 

<N ■ 

The probabilities of the two-photon spontaneous decay in the hydrogen atoms and hydrogenlike ions were studied since 
the theoretical formalism was worked out by Goeppert-Mayer [l| and the first estimate for the two-photon transitions 
O |' 2si/2 — > 2"/ (El) + lsi/2 was obtained by Breit and Teller |2|. This transition in H atom and low-Z H-like ions is of 
^ . special importance since it defines the lifetime of 2s level and exceeds by many orders of magnitude, approximately 
by (1/aZ) 4 , the probability of the one-photon decay 2s — > j(Ml) + Is. Here a is the fine structure constant, Z is the 
charge of the nucleus. A highly accurate nonrelativistic calculations of the transition probability for 2s — > 27(£1) + Is 
process for the hydrogen atom was performed by Klarsfeld 0. The first fully relativistic calculation of this transition 
probability was made by Johnson jj| and later by Goldman and Drake Q , , an d by Parpia and Johnson • The 
recoil corrections were given in the papers by Fried and Martin p| and by Bacher More recently, Karshenboim 
and Ivanov |Hj evaluated the radiative corrections to this decay. The decay probabilities ns — > 27(.ET) + Is with 
n = 3 — 6 in H atom were evaluated in [Til ] . 
' The E1M1 two-photon decay rate was so far evaluated only for the transition 2 3 Pq — > "/(El) + 7(M1) + l 1 5o in 
Q-ij two-electron ions. The reason is that for the level 2 3 Pn this decay channel is dominant in the absence of the hyperfine 
quenching as was first stated in [l2T ]. According to |l2j |. the angular momentum coupling rules for — > transitions 
t—{ | allow only the emission of two photons with equal values of the angular momentum. The probability of two-photon 
J> ■ decay 2 3 Po - * l(El) + 7(M1) + l 1 5o was evaluated within fully relativistic approach for Z = 92 in 0] an d f° r 
50 < Z < 94 in [14(. Recently a rigorous QED approach was applied to the evaluation of E1M1 transition in 
He-like ions with 30 < Z < 92 [Tr|. 

In this paper we describe the QED theory of two-photon decay process and calculate the E1M1 and E1E2 transition 
probabilities for H-like ions for nuclear charge Z values within the region 1 < Z < 100. 

Unlike the case of two-electron ions, both E1M1 and E1E2 transitions are allowed in the one-electron systems. As far 
as we know the transitions 2p — > j(El) +-f(Ml) + Is and 2p — > 7(-El) +7(_E2) + Is were never calculated. The reason 
is, of course, that they are very small (about 10 -13 for Z — 1) compared to the leading transition 2p — > 7(£T) + Is. 
Still both two-photon transitions behave like Z 8 for larger Z values. The Z behavior of El transition is essentially 
"Jj ■ weaker (Z ) and for Z = 92 the El transition prevails only by two orders of magnitude. 

We performed the calculation in two different gauges, thus receiving an accurate check of the gauge invariance of the 
, results. The two gauges employed were relativistic counterparts of the "velocity" and "length" forms of the transition 
Oh 1 amplitudes in the nonrelativistic case. In the low-Z region in the "velocity" gauge the intermediate negative-energy 
states give nearly dominant contributions both to E1M1 and E1E2 transition probabilities. Contrary to this, the 

■ negative-energy states contribution is fully negligible in the "length" gauge for the low-Z values. Similar conclusions 
on the importance of the negative-energy states in "velocity" gauge for the calculations of the one-photon transition 
probabilities in neutral atoms with one valence electron within the Relativistic Many Body Perturbation Theory 
approach were made earlier by Savukov and Johnson |17| . 

In this paper we derive also an explicit expression for the negative-energy contribution to the decay probabilities 
E1M1 and E1E2 for the low-Z H-like ions in "velocity" gauge. 

In the QED calculation of the two-photon transition probabilities one needs to generate the complete Dirac spectrum. 
In this paper we used the Dirac-Coulomb wave functions. For the summation over complete Dirac spectrum the 
different powerful numerical methods were developed in the last decades: the finite basis set method 0], the i?-spline 
method an d the space discretization method [2(J. In this work we used the version of the B-spline method 
presented in j^. The relativistic units are used throughout this paper. 



2. QED theory of two-photon transitions. 
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The two-photon decay process A — » A' + 2j for the noninteracting electrons is represented by the Feynman graphs 
Fig.l. In this section we characterize the photons by the momentum k and the polarization e. According to the 

Feynman correspondence rules the S-matrix element S A , a A is equal to [23 , H3- 



S { A'l = H) 2 e 2 / d i x l d A x 2 (^ A ,(x 1 )^ 1 A; i (x 1 )S(x 1 x 2 ) 1 ^ 2 A; 2 (x 2 )^ A (x 2 )) 



(1) 



Here 



*„(xi)*„(x 2 ) 



iO) +uji 

— oo 

is the electron propagator where the sum runs over the Dirac spectrum for the electron in the field of the nucleus 



(2) 



* A {x) = ^ A (v)e- lEAt 
is the electron wave function, E A is the electron energy, 



(3) 



^L e (X) e i(kr ,,/) 



(4) 



is the wave function of the photon characterized by the momentum k and polarization vector (/i, A = 1,2,3,4), 
x = (r, i). For the real transverse photons 



A e,k( r ) e " 



-ILUt 



(5) 



Inserting Eqs. ©-© m Eq. QJ, integrating over time and frequency variables and introducing the amplitude Ua>a as 



Saa 1 = —2TriS(E A > + w + J — E a )Ua>a 



we obtain 



(6) 



U A'A 



27re 2 ^ (" A e,k)^'«( aA e',k')" A 



E n -E A + u)' 



where e = y/a is the electron's charge. 
Defining the transition probability as 



aW * a = 2ttS(E a -E a ,-u)- J) \U% a + U\ b , 



dk dk' 

(2^)3 (2^f 



and integrating over u> we obtain finally 



dW A ' A {u>\v,v',e,e') = e 



w'{E A -E' A ~J) 



E 

n 

E 



(aA* k ) A - n (aA*, k ,) nA 

E n -Ea + lu' 

( aA e',k'^'"( aA e,k)'^ 



E n - E A + w 



(27T) 



dvdu'duj' 



(7) 



(8) 



(9) 



where v = k/u> 
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It is more convenient to come over to the photon's wave functions characterized by the angular momentum and 
parity. For this purpose we will use the expansion of the linearly polarized wave in spherical harmonics p3 | , |23| : 

ee ,(kr )= £ A W (r) (eYW (k)) (10) 

J,M,\=0,1 

where Ajj f (r) is the electromagnetic vector potential, Yj^-(k) is the vector spherical function of the magnetic (A = 0) 
or electric (A = 1) type. The electric and magnetic vector potentials are: 



A«(r)=, J+1 



| \f 2 /+ 1 gJ+i( fcr ) eY JJ+iM( r ) - sj 2J + +i 9J - 1 ( fcr ) eY J J-im ( r ) + 
Gj (\l 2 !J + 1 gJ+1 ( fcr ) eY ^+ 1M ( r ) + \[ 2 ^^j9J-i(kr)eY*j J _ 1M (r) + igj^ 1 (kr)Y^ M 



(11) 



A^(r)=z J ffJ (fcr)eY* 7JjU (r). (12) 

Here g^cor) — (2tt) 3 / 2 —== Ji+i/2(wr), J„ is the Bessel function, YjxAf(k) arc the vector spherical functions, J, M 
are the photons angular momentum and its projection, Gj is gauge parameter defining gauge for the electromagnetic 

potentials. In our calculations we employ the "velocity" gauge (Gj = 0) and the "length" gauge (Gj = )• 

Note, that Eq. © corresponds to Gj = 0. 

After substitution of Eqs. I|10I12|I in Eq. @ we can perform the summation over polarizations and integration over 
the photon's angles using the formula 

J2 J dv(e*Y L JM (v))(eY L j* M (v)) = J du{u x Y L JM {y)){y x Y L j',* M ,(u)) = &jj.& ll .& mm > (13) 
Then the expression for the two-photon transition probability looks like: 

(Ma'a{u) = e 4 — — - f -duj'x 

[Zn ) 



E E 

X,M,J X',M',J' 



( A )* \ I \ C^)* \ 

\ OL - Pi -J l M l Lj')n A ij A il A im A i,nj n l n m n ( aA JMLj)nj n l n m n ,n A j A lAmA 

E„,i„ — E. 



nj n - UnAjA + U (14) 



_j_ \ aA -JM^) n A'3A' l A' m A'' n Jn l n m ri V a A ./' M' u' ) " jn In ™.i , 71 A j A I A 771 A 

Enj„ - E, lAJA + lo' 

Here we have replaced each electron subscript A by the standard set of quantum numbers uaJaIaitia, where n is 
the principal quantum number, j, m are the total electron angular momentum and its projection and I defines the 
parity of the state. 

In this work we calculate the total rate 

"«?-ir"^ <i5) 

where u> max =E A - E A > ■ 

3. The angular reduction 

The angular integration in the matrix elements in Eq. I|14|) can be performed in a standard way: 

( a A-JM )n A ,] A ,-l A ,m A ,,n A „l A „] A „m A „ = Sj A ,m A , j A i,m A „ , J,A{Cj A ,l A , ,j A „l A „ , J^n A , j A , ,n A „ j A „ , J (1*0 

where 



c _ f_i \-m A ,-M-l/1 ( 3 A' j A" J 

Z] A ,m Al , 3A „m A ,,,J,M — { t) I ^ fhAn M J , 



the symbol to denotes —to, and 



C j Al l Al ,j A »l A »,J,M 



V4tt(2J+1) 



LU' , JA" 



(1/2 -1/2) n(A) (^''^"' J )' 



n (A) (Z A , , Z A /, , J) = for odd (I A ' + I A" + J + A) 
n (A) (Z A <, U", ■/) = 1 for even (l A , + l A » + J + A), 



\j,f] = vW + W + T) 
For the radial integrals we use the notations similar to ones in |l3j: 



d(0) ( \ 

n A ,j A ,,n A „j A „,J\ 1 



u> 2J+ 1 



2tt ^/J(J + 1) 



+k A ,,)I+ 



R nl>jA>,n A „j A »,J^ 



/7TT{( fc ' 4 '-^") / m + ( J + 1 ) / / + i} 



J + l 



{(Au'-fcA"K/-i- J/ 7-i} + 



G 7 ((2 J + + (k A , - k A „)(l} +1 + 7+_r) - JJJ_! + (J + l)/,7 +1 ) 



if = J r 2 drgj(uir)(g A 'fA" ± fA'9A" 



F.i = J r 2 drgj(ajr)(g A 'gA" + I A' I A") 
Here g n j(r), f n j{ r ) are the upper and lower radial components of the Dirac wave function and 

k-l 1 

+ if j = l+l 

The total decay rate should be summed over the magnetic quantum number m A i and averaged over the 
quantum number m A . Then 



dW A >A(<S) = 



2w(2j A + 1) 



E E E 

X,J,X',J' MM' m A /m A 



dui'x 

■^n A tj A tl A i,nj n l n ,J' ( U ')Sj A im A ,d n rn n J'MiT^ ln ^ nA j AlA j(u})Sj nmn j AmAt jM 



njnln 



+(u <-► cj', J <r+ J', A <-> A' , M <-> M') 
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where 



T 



(A) 



n A jAlA,n B jBlB,J\ I ~ 3AlA,jBlB,J n A j A ,,n B j B ,J\ I 

The summation over the M, M' , rriA' , rriA, m n can be carried out using the sum rules for 3j-symbols: 



>(A) 



(27) 



J A' m A' ,JAmA>J' M 1 ^ 'j A t ra A / ,j n t m n J ,J M va n i ,j A ra A ,J'M' 



' M' Sj ,,m A i, j„ i m„ i , JM Sj , m , 



all projections 



2i„ + 1 



(28) 



^ ^ ^j^/ m A i ,J n Tn n ,JM Sj n rn n ,j A m A ,J' M' Sj A , m A i ,j n rm n > , J' M' ^j n t 



3 n > m n' >lAm A 



all projections 



JM - M) \Ja- J in 



(29) 

Inserting the expressions (|16JI - H22|I in Eq. Ij26|l and performing summations over all the projection indices we obtain 
finally the formula for dW a< A expressed through the various radial integrals: 



dW A 2 , 1 }=e A — -i Y (dW^ + dW^+dW^)dw' 

2tt{2ja + 1) . ' _. 



A, A', J, J' 



(30) 



where 



47r(2j„ 



1 y 



n(A) 

'"■A'jA'lA' ,nj n ln,JV 1 ^ nj n l n ,n A j aI A ,J' 



(A') 



(31) 



dW (2) = dW (1) (A <-> A', J <-> J', w w') 



(32) 



8tt 2 [ JA' J 3n> J 



E T n A \j A d A ,,n) n l n ,A W ) T nj]l. n ,,n A 3 A l A ,J' T nJj A ,l A ,,nj,J n ,J< ( W ') T ijll n ,n AJA l A , J M 

4. E\E\ transition probabilities. 

The QED results of the calculations of the E1E1 two-photon transition probabilities for H-like ions with nuclear 
charges Z = 1 . . . 100 in comparison with results from la are given in Table 1. For the summation over the entire 
Dirac spectrum in Eqs. (|31|I - H33I) the B-spline approach |2lj was applied. The number of the grid points was N = 50; 
the order of splines k — 9. The radial integration was performed by the Gauss method with 10 integration points. 
Changing the number of the grid points, the order of splines and the number of integration points, we estimate our 
numerical inaccuracy as 10 -3 . We solved the Dirac equations with the Fermi model for the charge distribution p{r) 
inside the nucleus: 

p(r) = P \ — - (34) 

v ; exp[{r-c)/a] + 1 ' 

with a = 0.5350fm, po is defined from the normalization condition and c deduced via the relation An p(r)r 4 dr = 
(r 2 ). Here (r 2 ) 1 / 2 is the root-mean-square nuclear radius. 

The results of our calculation agree well with the results of calculation in p| . In the nonrelativistic limit the results 
excibit the behavior 

Wffif 1 = a£ 1 E 8 \aZ)°a.u. (35) 
with af^ 1 = 0.001317. This result coincides with the accurate nonrelativistic value Q. 
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5. E1M1 transition probabilities. 

The numerical results for E1M1 2p 1 / 2 — lsi/2 transition probabilities are given both in the "velocity" and "length" 
gauges in Table 2. The details of the calculations are the same as for the E1E1 transitions. The contributions of 
the positive-energy, negative-energy and the "total" contributions are shown separately. One has to remember, that 
the "total" contributions includes also the interference contribution between the positive-energy and negative-energy 
states. 

According to the coupling rules for E1M1 2pi/ 2 — * lsi/2 decay the sets of the intermediate states in the sum in 
Eg. 1141 are nsi/ 2 , npi/2, n P3/2 an d nd^/ 2 - Of them, the states with n = 1,2 give the dominant contribution for small 
Z values in the "length" gauge. This contribution scales like 

,j,E1M1(+) /, w ElMl(total) /, ,,>, E1M1/ r,\8 ; OI! \ 

W 2pls (length) = W 2 ls (length) = a 2 ls (aZ)°a.u. (36) 



with afflj 1 = 2.907 • 1CT 5 (37) 

The negative-energy and the interference contributions are quite negligible for small Z values in the "length" gauge. 
The scaling law Eg. 1361 can be understood from the estimate (in relativistic units): 



W. 



£1M1(+) 



2pls 



(length) 0;^ 

7T 



E 

n(+) 



n)(n|| M ||f) 



AE n 



i|H|n)(n||d||f) ) 
AE ni J 



(38) 



where u>if is the transition frequency 2p — ls, AE n i ar the effective energy denominators, d and fi are the electric 
and magnetic dipole transition operators, (a\\ .. . \ \b) are the reduced, i.e. integrated over the angles matrix elements 
and the summation is extended over the Schrodinger equation solutions. The matrix elements for the electric dipole 
operator d = y/ar (r is the electron radius-vector in an atom) are of the order (z||d||n) ~ y/a/maZ r.u. The matrix 
elements of the magnetic dipole operator fx = y/as/m (s is the electron spin) (n||/x||/) are the order y/a/m if the 
principal quantum numbers for the (n\ and |/) states coincide. Otherwise they are zero in the nonrelativistic limit 
due to the orthogonality of the radial wave functions. Then, assuming u>if ~ m(aZ) 2 r.u., AE n i ~ m{aZ) 2 r.u., and 
taking into account the relation a 2 r.u. = 1 a.u. for the energy units, we arrive at the result Eq. I|36|) . 

The picture looks quite different in the "velocity" gauge. In this case the contributions of the intermediate states 
with n — Is (first term in Eq. (J3EJ) and n — %P (second term in Eq. (jSHJ) cancel fully with the contribution of the 
interference term. This cancellation was checked numerically within the accuracy of our numerical procedure. The 
remaining positive-energy matrix elements of the magnetic dipole operator are nonzero only with the introduction 



of the spin-orbit interaction, i.e. are of the order (n||/i||/) 
positive-energy remainder is of the order 



(aZ) 

m 



a.u. Then the total contribution of the 



W 2 pif (velocity) w (aZ) 12 a.u. (39) 

This seems to be in agreement with the estimate, given by Drake [l3( for the E1M1 transition in two-electron ions 
in the high Z region (Z > 27) where the influence of the interelectron interaction becomes less significant. However, 
we have to remind that Eq. (|39l) gives only the minor contribution to the total two-photon E1M1 2p — ls decay 
rate for small Z values in the "velocity" gauge. The major contribution arises in this case from the negative-energy 
intermediate states and scales like 



^If^^ (velocity) = W^ 1(t ° tal) (velocity) = a^V^a.u. (40) 

with the same a^ptf 1 coefficients as in Eq.(j2Zjl. The scaling behavior in Ea. lj4Q(l will be demonstrated explicitly in 
Section 7. 

6. E1E2 transition probabilities. 

The E1E2 transition probabilities by the order of magnitude are comparable with the E1M1 transition probabilities 
for all Z values (see Table 3). This transition probability was also evaluated in the "velocity" and "length" gauge 
within the same numerical approach. 
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For the small Z values the E E1E2 transition probability scales with Z in the same way, as W E1M1 : 

W 2 E p \ E2 =a E p \f(aZfa.u. (41) 

with a E p 1E 2 = 1.986 • 1CT 5 (42) 

In the "length" gauge exclusively the positive-energy intermediate states contribute to the result Eq. I|41l) . The 
scaling law for this contribution follows from the estimate similar to Eq. I|.'j8[l: 



W E1E2(+) O 2 9 



<i||Q 20 ||n>H|d||/) , (t\\d\\n}(n\\Q 20 \\f) 



^ 1 SE„i AE„ 



(43) 



4?ra 2 



where Q 2 o is the spherical component of the quadrupole electric transition operator Q 2m — y - r l2ro(^)- Here 

12m is the spherical function, dependent on the electron angular variables. The matrix elements of Q in Ea. (|43|) are 
of the order (i||Q||n) ~ ^Ja/m(aZ) 2 . Inserting this estimate in Ea. H43l) we arrive at the result Ea. (|41|) . 

In the "velocity" gauge all the contribution from the positive-energy, negative-energy and intermediate parts are 
comparable. The value and the scaling behavior of the negative-energy part will be derived analytically in Section 7. 

7. Analytic derivation of the negative-energy contribution to the E1M1 and E1E2 2pi/ 2 — lsi/2 
transition probabilities in the "velocity" gauge for small Z values. 

In this section we derive an explicit expression for the negative-energy contribution to the E1M1 and E1E2 2pl/2— lsj/2 
transition probabilities in the "velocity" gauge for small Z values. This derivation will help us to check the validity 
of our numerical calculations in the nonrelativistic region (small Z values). 

We will perform this derivation using another set of the photon's characteristics, namely photon's momentum k 
and polarization vector e. Thus we will not be able to distinguish between E1M1 and E1E2 transition probabilities 
and will evaluate them as a unique correction to the dominant El 2pi/ 2 — lsi/2 transition. 

The starting point for our calculations is the formula 10 where we retain the sum only over the negative-energy 
states. The corresponding energy denominators in the nonrelativistic regime we replace by -2m, neglecting also the 
photon frequencies, limited by the value u>if = Ea — Ea 1 — E 2pi/2 — E 2si/2 . Here m is the electron mass. We also 
expand both exponents Eq.Q in Eq.©, replacing one of these exponents by the unity and retaining omly the next 
term of the expansion in another exponent. After the summation over photon's polarizations and the integration 
over the photon's emission directions this will give us the desired correction to the leading El 2p 1 / 2 — lsi/2 transition 
amplitude. 



Thus we start with the expression 



(44) 



where 



U[j ] = ~ {(i|(e -)(kr)|n ( _ ) )(n ( _ ) |(e' CT )|/) + (i|(e < T)|n ( _ ) )(n ( _ ) |(e' < r)(k'r)|n ( _ ) ) + [e,k <-> e',k']} (45) 

Ea. l|45() corresponds to the "velocity" gauge. Here cr are the Pauli matrices and the summation is extended over the 
negative-energy set of solutions of the Dirac equation for the electron in the field of the nucleus. In the nonrelativistic 
limit this set comes over to the complete set of the solutions of the Schrodinger equation for the positron in the field 
of the nucleus. 

The employment of the closure relation then yields 



^_ ) = _ ! (eeO 
J m 



For the summation over the polarizations we apply the standard formula [2/ 



e 



(47) 
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Then 

^(ee')*(ee') = 2-(^') 2 (48) 
ee' 

For the integration over the photon's emission angles the following relations can be used: 

duvi = I dvviv^vi = (49) 



Air 

duuiv k = -g-^ife (50) 
Ait 

dvviVjVkVi = — (SijSki + 5i k Sji + SuSjk) (51) 
The summation over the polarization and the integration over the photon's emission angles in Eg . (|44|l yields 

dvdv'(ki + kl)(kj + fc')[2 - (y x i/') 2 ] = {uj 2 + uj' 2 )—^-5 l] (52) 

9 

Then, inserting Ea. (|o"2|> in Ea.l|44|) and using Eq. l(T5|) we obtain 

K~) _ " l/„-|».IA|2 / , ,/, , , ,U, ,2 i (.. , ,2\2i j, , _ zu , ,5 f\|2 



W « = 9^ mf)l J ^f-^ + ("if-» 2 r]d» = ^v^/K^I/)! 2 (53) 
The order of magnitude and the scaling behavior for the transitions 

<!f (velocity) + W^f/" (velocity) = ^ (aZ)\.u. (54) 



with 



a[ p \ n s r = 5.625 • 10~ 5 (55) 



follow immediately from Ea. (|53|l in the nonrelativistic limit. This result can be compared with the exact result of the 
numerical evaluations: 

Ku (velocity) + Wg£ 2( -> (velocity) = a$ B (aZ) 8 a.u. (56) 

with 



a 



i~l = 5.806 • 10- 5 (57) 



The discrepancy between a^piL &nd o^pis (3,1%) is larger than the expected discrepancy due to the relativistic 
corrections 



(-) „(-)™-\ 

{aZf (58) 



(a\-> -a ( ~ ,nr ) 

\ a 2pls u 2pls ) 



a 2pls 

for Z=l. Besides, the discrepancy in our case does not depend on Z for small Z values. It can be argued, however, that 
the order of magnitude and Z-dependence l|58|) correspond to the positive-energy contributions and can be different 
in case of the negative-energy intermediate states. 

In total, the discrepancy between a 2p [ s and a 2p [ s is small enough to confirm the validity of our numerical calcu- 
lations in the region of the small Z values. 

8. Conclusions. 
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In this work we have evaluated for the first time the two-photon transitions probabilities W^^ 1 and W^\f 2 in the 
H atom and in H-like ions with Z up to Z = 100. The evaluation is performed in a fully relativistic way. For the 
small Z values the scaling law (aZ) 8 is established for both probabilities. The total transition probability for 2p — > Is 
transition can be presented in a form 



W = W, 



7T 



1 + ~{aZ) 4 F(aZ) a.u. (59) 



where Wq is the transition probability for 2p — > "/(El) + Is process: 

W Q = ^f f \(t\r\f)\ 2 a 3 a.u. (60) 

The function F(aZ) for all Z is given in Table 4. For Z = 1 the correction term in Ea. i|59|) is quite small, but grows 
very fast with the increase of the nuclear charge Z and becomes significant in the one-electron highly charged ions. 
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Table 1. Total two-photon E1E1 2si/ 2 — lsi/2 transition probabilities in sec 1 for different Z values. Numbers in 
parentheses are the powers of 10. 
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yy \ v '"V 


TT/"('lcn) 
yy v iL,i v 
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yy \ vv -v 


yy v iv "-v 


TT/" [5l 
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O 1 A CT O { 1 A\ 
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tr o/ 1 A t: /oA 
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c o/^at/oN 

5.2607(2) 
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A C 
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O A AAO / 1 A\ 

6.4003(10) 
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c aaaa/o\ 

5.9909(d) 


5.9910(3) 


r: aato /o\ 
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TZ A 

50 
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1 1 OC /I /1 1 \ 

1.1854(11) 


1 1 OZ? A / 1 1 \ 

1.1869(11) 


4 


3.3652(4) 


3.3653(4) 


3.3689(4) 


56 


O OA /i *7 ( 1 1 \ 

2.2947(11J 


O OA AO 1 A 1 A 

2.2948(11) 


O OAO A / 1 1 \ 

2.2980(11) 


5 


1.2833(5) 


1.2834(5) 


1 no ^/r\ 

1.2847(5) 


60 


O 1 O O A /' 1 1 \ 

3.4229(11) 


3.4230(11) 


3.4282(H) 


6 


3.8305(5) 


3.8306(5) 


3.8347(5) 


65 


5.4293(H) 


5.4293(11) 


5.4387(11) 


7 


a / i rrn/r\ 

9.6550(5) 


9.6551(5) 


9.6654(5) 


70 


8.2975(11) 


8.2975(11) 


8.3139(11) 


8 


o 1 rnri// 1 ) 

2.1500(6) 


o i r A o / z" 1 \ 

2.1502(6) 


2.1525(6) 


74 


1.1379(12) 


1 1 O TA / 1 O \ 

1.1379(12) 


1.1404(12) 


Q 




4;.OtJU0 J 




80 


1.7655(12) 




1 77D1 f1 9 s ! 


10 


8.1921(6) 


8.1922(6) 


8.2010(6) 


85 


2.4747(12) 


2.4748(12) 


2.4824(12) 


12 


2.4425(7) 


2.4425(7) 


2.4451(7) 


90 


3.3899(12) 


3.3899(12) 


3.4021(12) 


15 


9.2914(7) 


9.2915(7) 


9.3017(7) 


92 


3.8216(12) 


3.8216(12) 


3.8361(12) 


20 


5.1898(8) 


5.1899(8) 


5.1956(8) 


100 5.9782(12) 


5.9783(12) 


6.0045(12) 


30 


5.8151(9) 


5.8152(9) 


5.8217(9) 
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Table 2. The two-photon E1M1 2p 1 / 2 — lsi/2 transition probabilities in sec 1 for different Z values in the "velocity" 
and "length" gauge. The positive-energy W + , negative-energy W_ and total contributions W t are given. 



V 

Zi 


velocity 


length 




W+ 


W- 


w t 


W+ 


W- 


w t 


1 


3.223(-5) 


9.667(-6) 


9.667(-6) 


9.667(-6) 


5.889(-16) 9.667(-6) 


2 


8.249(-3) 


2.475(-3) 


2.474(-3) 


2.474(-3) 


2.410(-12) 2.474(-3) 


4 


2.111 


6.334(-l) 


6.332(-l) 


6.331(-1) 


9.886(-9) 


6.331(-1) 


8 


5.407(2) 


1.625(2) 


1.619(2) 


1.619(2) 


4.056(-5) 


1.619(2) 


10 


3.219(3) 


9.702(2) 


9.637(2) 


9.637(2) 


5.910(-4) 


9.637(2) 


14 


4.745(4) 


1.437(4) 


1.418(4) 


1.418(4) 


3.361(-2) 


1.418(4) 


18 


3.538(5) 


1.078(5) 


1.054(5) 


1.051(5) 


6.888(-l) 


1.055(5) 


20 


8.211(5) 


2.510(5) 


2.443(5) 


2.433(5) 


2.445 


2.443(5) 


24 


3.523(6) 


1.086(6) 


1.045(6) 


1.039(6) 


2.193(1) 


1.045(6) 


28 


1.206(7) 


3.754(6) 


3.562(6) 


3.533(5) 


1.404(2) 


3.563(6) 


30 


2.092(7) 


6.547(6) 


6.164(6) 


6.105(6) 


3.227(2) 


6.164(6) 


34 


5.674(7) 


1.798(7) 


1.664(7) 


1.643(7) 


1.462(3) 


1.664(7) 


38 


1.376(8) 


4.421(7) 


4.014(7) 


3.951(7) 


5.609(3) 


4.014(7) 


40 


2.070(8) 


6.700(7) 


6.020(7) 


5.914(7) 


1.044(4) 


6.021(7) 


44 


4.417(8) 


1.453(8) 


1.276(8) 


1.249(8) 


3.315(4) 


1.277(8) 


48 


8.812(8) 


2.952(8) 


2.530(8) 


2.465(8) 


9.543(4) 


2.531(8) 


50 


1.218(9) 


4.119(8) 


3.486(8) 


3.386(8) 


1.569(5) 


3.486(8) 


54 


2.240(9) 


7.734(8) 


6.366(8) 


6.151(8) 


4.011(5) 


6.366(8) 


58 


3.939(9) 


1.391(9) 


1.111(9) 


1.068(9) 


9.614(5) 


1.112(9) 


64 


8.549(9) 


3.134(9) 


2.386(9) 


2.269(9) 


3.221(6) 


2.387(9) 


68 


1.375(10) 


5.180(9) 


3.813(9) 


3.597(9) 


6.807(6) 


3.813(9) 


70 


1.725(10) 


6.593(9) 


4.767(9) 


4.479(9) 


9.744(6) 


4.767(9) 


74 


2.659(10) 


1.048(10) 


7.307(9) 


6.805(9) 


1.943(7) 


7.308(9) 


78 


3.999(10) 


1.629(10) 


1.094(10) 


1.009(10) 


3.747(7) 


1.094(10) 


80 


4.863(10) 


2.016(10) 


1.327(10) 


1.218(10) 


5.145(7) 


1.328(10) 


84 


7.075(10) 


3.042(10) 


1.927(10) 


1.749(10) 


9.502(7) 


1.927(10) 


oo 


1.009(11) 


4.511(10) 


2.748(10) 


2.466(10) 


1.711(8) 


2.748(10) 


90 


1.196(11) 


5.460(10) 


3.262(10) 


2.910(10) 


2.276(8) 


3.262(10) 


92 


1.411(11) 


6.583(10) 


3.859(10) 


3.422(10) 


3.012(8) 


3.859(10) 


94 


1.658(11) 


7.910(10) 


4.550(10) 


4.010(10) 


3.967(8) 


4.551(10) 


98 


2.260(11) 


1.130(11) 


6.269(10) 


5.456(10) 


6.781(8) 


6.270(10) 


100 


2.621(11) 


1.344(11) 


7.329(10) 


6.339(10) 


8.805(8) 


7.330(10) 
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Table 3. The two-photon E1E2 2p 1 / 2 — lsi/2 transition probabilities in sec 1 for different Z values in the "velocity" 
and "length" gauge. The positive-energy W + , negative-energy and total contributions W t are given. 



z 


velocity 


length 




W+ 


W- 


w t 


W+ 


W- 


W t 


1 


1.232(-6) 


9.667(-6) 


6.605(-6) 


6.604(-6) 


3.625(-27) 6.604(-6) 


2 


3.154(-3) 


2.474(-3) 


1.690(-3) 


1.691(-3) 


2.375(-22) 1.691(-3) 


4 


8.075(-2) 


6.335(-l) 


4.327(-l) 


4.326(-l) 


1.557(-17) 4.326(-l) 


8 


2.065(1) 


1.620(2) 


1.105(2) 


1.106(2) 


1.023(-12) 


1.106(2) 


10 


1.230(2) 


9.655(2) 


6.583(2) 


6.582(2) 


3.638(-ll) 


6.582(2) 


14 


1.813(3) 


1.423(4) 


9.683(3) 


9.681(3) 


7.949(-9) 


9.681(3) 


18 


1.351(4) 


1.061(5) 


7.197(4) 


7.196(4) 


4.452(-7) 


7.196(4) 


20 


3.136(4) 


2.463(5) 


1.668(5) 


1.667(5) 


2.409(-6) 


1.667(5) 


24 


1.345(5) 


1.057(6) 


7.125(5) 


7.125(5) 


4.482(-5) 


7.125(5) 


28 


4.605(5) 


3.620(6) 


2.428(6) 


2.427(6) 


5.321(-4) 


2.427(6) 


30 


7.987(5) 


6.278(6) 


4.198(6) 


4.198(6) 


1.611(-3) 


4.198(6) 


34 


2.167(6) 


1.703(7) 


1.132(7) 


1.132(7) 


1.206(-2) 


1.132(7) 


38 


5.262(6) 


1.135(7) 


2.728(7) 


2.728(7) 


7.227(-2) 


2.727(7) 


40 


7.920(6) 


6.222(7) 


4.088(7) 


4.088(7) 


1.652(-1) 


4.088(7) 


A A 

44 


1.692(7) 


1.328(8) 


8.654(7) 


8.655(7) 


7.696(-l) 


8.654(7) 


48 


3.386(7) 


2.654(8) 


1.712(8) 


1.712(8) 


3.144 


1.712(8) 


50 


4.688(7) 


3.671(8) 


2,356(8) 


2.356(8) 


6.093 


2.355(8) 


54 


8.658(7) 


6.762(8) 


4.290(8) 


4.291(8) 


2.126(1) 


4.290(8) 


58 


1.530(8) 


1.191(9) 


7.465(8) 


7.468(8) 


6.807(1) 


7.465(8) 


64 


3.361(8) 


2.596(9) 


1.592(9) 


1.594(9) 


3.404(2) 


1.592(9) 


68 


5.460(8) 


4.189(9) 


2.530(9) 


2.532(9) 


9.217(2) 


2.530(9) 


70 


6.888(8) 


5.265(9) 


3.153(9) 


3.155(9) 


1.486(3) 


3.152(9) 


74 


1.077(9) 


8.152(9) 


4.793(9) 


4.801(9) 


3.726(3) 


4.793(9) 


78 


1.646(9) 


1.232(10) 


7.103(9) 


7.115(9) 


8.941(3) 


7.103(9) 


80 


2.020(9) 


1.502(10) 


8.568(9) 


8.589(9) 


1.364(4) 


8.569(9) 


84 


o.UUz^y ) 


z.iyy j 


i.ZZuyW ) 


l.zzof 1U ) 


3.092(4) 


1.226(10) 


oo 


4.390(9) 


3.158(10) 


1.717(10) 


1.723(10) 


6.781(4) 


1.717(10) 


90 


5.280(9) 


3.759(10) 


2.017(10) 


2.025(10) 


9.927(4) 


2.016(10) 


92 


6.329(9) 


4.455(10) 


2.357(10) 


2.367(10) 


1.444(5) 


2.356(10) 


94 


7.562(9) 


5.259(10) 


2.740(10) 


2.754(10) 


2.085(5) 


2.740(10) 


98 


1.070(10) 


7.239(10) 


3.653(10) 


3.677(10) 


4.272(5) 


3.653(10) 


100 


1.268(10) 


8.444(10) 


4.187(10) 


4.218(10) 


6.060(5) 


4.187(10) 



Tabic 4. The function F(aZ) ■ 1CT 3 for different Z values. 



Z F(aZ) 


Z 


F(aZ) 


1 


3.945 


50 


3.569 


2 


3.943 


54 


3.508 


4 


3.942 


58 


3.443 


8 


3.935 


64 


3.340 


10 


3.929 


68 


3.268 


14 


3.915 


70 


3.231 


18 


3.894 


74 


3.153 


20 


3.883 


78 


3.075 


24 


3.857 


80 


3.034 


28 


3.825 


84 


2.954 


30 


3.807 


88 


2.873 


34 


3.768 


90 


2.833 


38 


3.725 


92 


2.794 


40 


3.702 


94 


2.754 


44 


3.651 


98 


2.677 


48 


3.597 


100 


2.640 
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FIG. 1: The Feynman graphs, corresponding to the two-photon decay process A — > A' . The photons are characterized by the 
momentum k and the polarization e 



